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Consider a ring of iV qubits in a translationally invariant quantum state. 
We ask to what extent each pair of nearest neighbors can be entangled. 
Under certain assumptions about the form of the state, we find a formula 
Q . for the maximum possible nearest-neighbor entanglement. We then compare 

this maximum with the entanglement achieved by the ground state of an 
antiferromagnetic ring consisting of an even number of spin- 1/2 particles. 
We find that, though the antiferromagnetic ground state typically does not 



maximize the nearest-neighbor entanglement relative to all other states, it 



does so relative to other states having zero ^-component of spin. 



PACS numbers: 03.67.-a, 03.65.Bz, 75. 10. J 



1 Introduction: Entanglement Sharing 

Quantum entanglement, as exemplified by the singlet state of two spin- 1/2 
particles, tttjGTI) - lit)); has been the subject of much study in recent years 
[|], largely because of its connection with quantum communication @] and 
computation H. Entanglement bears some resemblance to classical correla- 
tion, but it differs in important respects, including the fact that entangled 
objects can violate Bell's inequality ||. Perhaps one of the most charac- 
teristic differences is this: if two similar quantum objects are completely 
entangled with each other, then neither of them can be at all entangled with 
any other object, whereas there is no such restriction on classical correlations. 
This property is sometimes called the "monogamy" of entanglement. For the 
special case of three binary quantum objects — three qubits — a quantitative 
extension of this rule has been proven in terms of a measure of entanglement 
called the "concurrence" which takes values between zero and one: the square 
of the concurrence between qubits A and B, plus the square of the concur- 
rence between qubits A and C, cannot exceed unity ||. In other words, to 
the extent that qubits A and B are entangled with each other, they limit the 
entanglement between qubits A and C. 

The present paper further explores the degree to which entanglement 
can be shared among a number of qubits. We focus on two closely related 
but distinct problems, (i) We consider a ring of N qubits in a translationally 
invariant pure quantum state and ask to what extent nearest neighbors can be 
entangled with each other; specifically, we ask how large the nearest-neighbor 
concurrence can be. Note that in this first problem there is no Hamiltonian 
specified; we are simply asking about the entanglement characteristics of 
quantum states, (ii) For our second problem we consider a particular physical 
system, namely a ring of N spin-1/2 particles interacting via the Heisenberg 
antiferromagnetic Hamiltonian, and ask whether the ground state of this 
system is a state of maximum nearest-neighbor entanglement. We will find 
that the antiferromagnetic ring maximizes entanglement within a limited set 
of states, but not absolutely. 

In both of these problems, we are focusing on pairwise entanglement 
within a system of iV particles. At least three problems with a similar focus 
have been considered before. Diir pi has shown that given a system of N 
qubits and any specified set of pairs of those qubits, one can design a state 
such that all the pairs in the chosen set are entangled and all the other 



pairs are not. Koashi et al. Jl]] have studied completely symmetric states of 
TV qubits and have found that the maximum possible concurrence between 
pairs is exactly 2/N. Thus in this context where all the qubits are required to 
be equally entangled with each other, the pairwise entanglement goes to zero 
in the limit of an infinite collection. Wootters [§] has considered a different 
problem, in which the qubits are arranged in an infinite line and only the 
nearest-neighbor entanglement is maximized. He found that for the infinite 
chain in a translationally invariant state, the nearest-neighbor concurrence 
does not have to be zero but can be as large as 0.434. It is not yet known 
whether this value is optimal. The problem we are about to address is the 
simplest finite version of the infinite chain problem. 

There have been several other studies of entanglement in iV-component 
systems, usually focusing on higher-order rather than pairwise entanglement 
M. All of these studies contribute to our understanding of entanglement 
distributed among many objects. We hope that our present results can 
eventually be combined with other work to construct a general theory of 
entanglement-sharing, not limited to qubits or to any particular geometry. 

2 Maximizing Nearest-Neighbor Entanglement 



We begin by recalling the definition of the concurrence |L0], |Tl|] between a pair 
of qubits, which we will think of as spin-1/2 particles. Let p be the density 
matrix of the pair, expressed in the standard basis {|TT)? |TI)> lit); III)}- Let 
p be the spin-reversed density matrix, defined by p = (a y ® a y )p T (a y <S> cr y ), 

where a y is the matrix . and the superscript T indicates transposi- 

tion. Then the concurrence of p is given by C = max{Ai — A 2 — A 3 — A4, 0}, 
where the Aj are the square roots of the eigenvalues of pp in descending order. 
(These eigenvalues are guaranteed to be real and non-negative even though 
pp is not necessarily Hermitian.) Concurrence is justified as a measure of en- 
tanglement by a theorem JTI| showing that C is a monotonically increasing 



function of the entanglement of formation [1^]. As we mentioned above, the 
values of concurrence range from zero, for an unentangled state, to one, for 
a completely entangled state such as the singlet state. 

We imagine a set of N particles arranged in a ring, with the locations of 
the particles labeled by an integer i — 1, . . . , N. In defining our problem, 



we restrict our attention to translationally invariant pure states \ijj), that is, 
states that under the cyclic permutation i — > i + k (mod N) are changed by 
at most an overall phase factor. This restriction forces the concurrence to be 
the same for each pair of nearest neighbors. The problem, then, is simply to 
find the maximum possible value of this concurrence .[] 

We have not yet been able to solve this problem in full. We solve instead 
a more tractable problem in which we limit the set of states over which the 
maximization is to be done. Specifically, we require our states to satisfy the 
following two conditions:^ 

1. The state \ip) of the ring is an eigenstate of the total z-component of 
spin. 

2. Neighboring particles cannot both be in the state ||). 

Though we are clearly leaving out many possible states, it is plausible that 
the maximum value we obtain for our restricted problem will not be far from 
the absolute maximum. This is because our two conditions tend to favor 
states with high nearest-neighbor entanglement. To see this, let us consider 
the density matrix p of a pair of nearest neighbors, obtained by tracing \ip) (i/j\ 
over all the other particles. Condition 1 implies that for any pair of particles, 
there can be no coherent superposition of basis states with different numbers 
of up-spins, e.g., |||) and |||), because the corresponding states of the rest 
of the chain are orthogonal. The density matrix p must therefore have the 

1 For a general, non-translationally-invariant state, one could define at least two distinct 
problems along similar lines: (i) maximize the average entanglement over all nearest- 
neighbor pairs, and (ii) maximize the minimum entanglement of all nearest-neighbor pairs. 
The first of these problems could be sensitive to the measure of entanglement one is 
using — e.g., concurrence, squared concurrence (also called the "tangle"), or entanglement 
of formation — even though these are all monotonic functions of each other. Problem (ii), 
which does not have this sensitivity, may thus be more interesting; it may also reduce to 
the translationally invariant problem considered here. 

2 Condition 2 breaks the symmetry between ||) and ||). Our choice to use ||) rather 
than ||) in the statement of this condition is arbitrary and does not affect any of our 
conclusions. 



following block diagonal form:[] 



f v \ 
w z 
0^x0 

V0 y J 



(1) 



One can show by direct calculation that the concurrence of this density ma- 
trix is 

C = 2msLx{\z\-^/vy,0}. (2) 

Condition 2 implies that the matrix element v is zero, so that the neighboring- 
pair density matrix becomes 



( ^ 

w z 

z x 

V0 yj 



(3) 



and the concurrence becomes simply 



C = 2\z\ 



(4) 



Density matrices of the form (|3|) have been singled out in two recent stud- 
ies as having particularly high entanglement. Specifically, Ishizaka and Hi- 
roshima [|n| have proven that such density matrices maximize entanglement 
for a fixed set of eigenvalues when one of the eigenvalues is zero. (They also 
show numerically that the form ([!]) is optimal when all four eigenvalues are 
non-zero.) Munro et al. |TJJ have shown that certain states of the form 
(||) maximize concurrence for a fixed value of the purity, defined as Tr(p 2 ). 
These studies suggest that our two conditions are consistent with high entan- 
glement, but they do not guarantee that we will be able to reach the absolute 
maximum. Indeed, we will see below that for at least one value of N, the 
optimal concurrence is not achievable by any state satisfying our conditions. 
Nevertheless, our solution to the restricted problem will be useful in Section 3 



3 In fact translational invariance implies that the matrix elements w and x must be 
equal — the frequency of occurrence of |||) in the ring must be the same as that of |||) — 
but we will not need to use this equality in what follows. 



where we discuss antiferromagnetic rings, and it should also serve as a good 
starting point for future work on the complete problem. 

Condition 1 forces the ring's state \ip) to have a fixed number p of up- 
spins and a fixed number N — p of down spins, but it does not specify the 
value of p. Our strategy will be to fix the values of both N and p and to 
maximize the nearest-neighbor concurrence over all states having these values 
and satisfying condition 2. This problem turns out to be exactly soluble, so 
that one can write down an analytic formula for the maximum concurrence 
Cmax{N,p). We can then use this formula to find the optimal number of 
up-spins, and hence the optimal concurrence, for any ring size N. 

For fixed N and p, the most general state we are considering has the form 

l^)= S K...i p \h---ip), (5) 

l<il<—<i v <N 

where \i\. . . i p ) is the state in which the particles at locations i\, . . . , i p have 
their spins up and all the other particles have their spins down. Though the 
above sum requires values of b only for sets of indices that are in ascending 
order, for convenience we define b to be symmetric in all its indices and equal 
to zero if any two indices have the same value. The normalization condition 
on the coefficients b is 

E IW = i- ( 6 ) 

l<il<-<ip<JV 

The condition of translational invariance is expressed as 

\ — p ikd h 

H\...i v C "i 1 +fc...ip+fc; 



h x ...i p — e 4 bi 1+ k...ij,+k, (7) 



where addition is understood to be mod N and e lNe = 1. Finally, in accor- 
dance with condition 2 above, the coefficients b must satisfy the constraint 

bi x ...i p = if i n -i m = l for any n,m= 1,. . .,p. (8) 

That is, no state is allowed in which two up-spins are adjacent. 

To find the concurrence between two neighboring particles, we need to 
find the off-diagonal element z of the two-particle density matrix as expressed 
in Eq. @. Translational invariance guarantees that the value of z will be 
the same for each pair of nearest neighbors; we consider a specific pair at 
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locations i and % + 1. Taking the partial trace of \ip){ip\ over all the other 
particles, we find that 

z = /2 h,k2...k p bi+i,k 2 ...k p , (9) 

l<k 2 <-<k p <N 



so that 



C = 2\z\ 



/2 2bito...k p bi+i,k 2 ...k p ■ (10) 

l<fc 2 <-<fc p <Af 

This form tells us immediately that the concurrence can be maximized by 
choosing the coefficients b to be real and non-negative: if we were to use 
complex values, then the concurrence could only be made larger, not smaller, 
by replacing each coefficient b by its absolute value. Let us therefore restrict 
our attention to such real and non-negative states. In that case, translational 
invariance takes the simple form 

Ou..i p = Ou+fc,...,i p +fc- (11) 

Thus, once the values of fri,^..., are fixed, all the other fe's are determined. 

The condition expressed by Eq. (|8|), i.e., that no two up-spins should be 
adjacent, is an awkward one to enforce directly. It is therefore helpful to relate 
our problem to a simpler problem that does not have this constraint. Roughly 
speaking, we do this by removing from the ring the site immediately to the 
right of each up-spin. More precisely, we consider a ring of N — p particles 
with exactly p up-spins, and we assign to every state \tp) of our original ring 
(every state, that is, that satisfies our conditions) a corresponding state \(p) 
of the smaller ring: 

10)= E d n .., p \ 3l ... Jp ). (12) 

l<ji<-<j p <N-p 

The coefficients d are defined in terms of the original coefficients b\^,„i v with 
1 < % 2 < ■ ■ ■ < i p . Let j 2 = h - 1, J3 = »3 - 2 , • • • ,J P = i P ~ ip ~ 1); then 
dxj 2 ,„j p = JN/(N — p)b\ t i 2 .„i p - The values of the other cfs are determined 
by translational invariance — that is, dj x .„j v = dji+fc,...j„+ifc (mod iV — p)— 
and as before, we take dj lm * to be symmetric under permutations of the 
indices and equal to zero whenever two indices have the same value. The 



factor JN/(N — p) is included in order to make \<f>) normalized: translations 
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around the ring generate fewer efs than 6's, so that the cfs need to be larger .0 
Let us define a pseudo-concurrence C of the smaller ring by analogy with 
Eq. (0). 

C = 2^ 2dj t k 2 ...k p dj+i > k 2 ...k p , (13) 

l<fc2<-<fcj»<iV-p 

where we have omitted the absolute value sign since the cfs are all real and 
non-negative. Because our states of the small ring do not satisfy condition 2, 
C is not the nearest-neighbor concurrence of the state \<f>). However, because 
of the relationship between d and b, we can use C to find the concurrence C 
of our original ring: 

C=^C. (14) 

Thus we want to find the maximum possible value of C over all real and 
non-negative, translationally invariant states of the (N — p)-particle ring 
with exactly p up-spins. 

To do this, let us rewrite Eq. ([HJ) in a more convenient form by intro- 
ducing the creation and annihilation operators 

a ) = f o o ) and aj = ( i o J ' (15) 

which act on the particle at location j of the small ring and are expressed 
here in the basis {|T)> 14)}- I n terms of these operators, Eq. (|T3D becomes 
simply 

C' = 2{<P\a) +1 a^). (16) 

The value given by Eq. (p!6|) is the same for all pairs {j,j + 1}. We can 
therefore write C as the average of this quantity over j: 

9 N-p 

C' = -—{<j>\J2aUa j \(f>). (17) 

iv p j=1 

Again using our assumption that the coefficients are real, we can re-express 
Eq. (0) as 

c ' = -(^)»). ( 18 ) 

4 For each collection of 6's that are equal because of translational invariance, there is a 
corresponding set of d's, and the ratio of the sizes of these sets is always N/(N — p). 



(19) 



where 

N-p 

H = -^2 (4 a i+i + a )+\ a j)- 
i=i 

In other words, a state |0) maximizes C if it minimizes the expectation value 
of the operator H, as long as this minimum is achieved with only non- negative 
real values of the coefficients d. 

The operator H is the Hamiltonian for the one-dimensional ferromagnetic 
XY model; so our problem reduces to finding the lowest-energy state of 
this model with exactly p spins up. This is a solved problem JT5[]. The 
solution begins with the observation that the operators at and a are not 
quite fermionic creation and annihilation operators, since [a,-, Ofe] = [%, a k ] = 
[a], a k ] = for j ^ k, whereas truly fermionic operators attached to different 
sites would anticommute. It is helpful to define new creation and annihilation 
operators $ and c that are genuinely fermionic: 



j'-i 



cxp 



[ i7r E 



a\a k 



k=\ 



"3i 



(20) 



J 



olexp 



j-i 



in Y^ a k a k 
fc=i 



In terms of the c operators, we have 

N-p 

J. 



H — — y^ (cjcj + i + cJ +1 Cj) if p is odd 



(21) 



(22) 



and 



(JV-p)-l 

H = - Yl ( c ] c j+i + c ]+i c j) + ( C N- P c i + c\cn- p ) if P is even. 



(23) 



For either odd or even p, the Hamiltonian can be diagonalized exactly, so 
that the system can be regarded as a collection of p independent identical 
fermions. For odd p, one finds that the energy eigenvalues of these fermions 



are e r . 



-2 cos 



—2 cos 



( (2to+1> 



N-p 



m = 1, . . . , N — p, whereas for even p they are 

y v _ i; I, m = 1, . . . , iV — p. The minimum value of ((p\H\(p) is 
the sum of the p smallest values e m , since in the ground state the fermions 



will occupy the p lowest energy levels. This sum turns out to be given by the 
following formula, valid for both even and odd values of p. 

2sin( 1 |^ N ) 
E ■ - KN ~ pJ (241 

sm (iv^J 

The state \<f>) corresponding to this energy is the discrete version of the 
ground-state wavefunction of a set of p hard beads on a loop of wire. The 
coefficients dj 1 ...j associated with this state can be taken to be real and non- 
negative, and the state is translationally invariant. Thus the assumed condi- 
tions are met and we can use E min to find the maximum pseudo-concurrence 
C max in accordance with Eq. (pi 



1 2sin(^) 

N-p (Ar-^sm^; 



Finally, using the relation (|1J), we get the maximum nearest-neighbor con- 
currence of our original ring of N particles: 

2 sin (^-) 

C max (N,p) = y^-. (26) 

Again, this is the maximum value under the following assumptions: (i) the 
ring has exactly p spins up, and (ii) no two up-spins are adjacent. 

For a given value of N, we now need to find out what value of p maximizes 
C m ax{N,p). For any fixed N it is easy enough to carry out this maximization 
explicitly. Consider, for example, the case N = 7. In a ring of seven parti- 
cles, the number p of up-spins can have any of the following values without 
violating our condition 2: p = 0,1,2, and 3. Inserting these numbers into 
Eq. ( p6|) we get the corresponding values of the concurrence: C = 0, 0.286, 
0.462, and 0.286. Thus for a ring of seven particles it is best (under our 
assumptions) to have two spins up and five spins down. We have carried out 
this sort of direct maximization for the first several values of the ring size N, 
with the following results: 



N 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Popt 


1 


1 


1 


1 or 2 


2 


2 


2 


3 


3 


^max 


1.000 


0.667 


0.500 


0.400 


0.471 


0.462 


0.433 


0.444 


0.449 



10 



Note that though the maximum concurrence tends to decrease with increas- 
ing N, it is by no means monotonic. 

It is interesting to find the limiting value of C max as iV goes to infinity. 
To do this, we write Eq. ([26|) in terms of iV and a = p/N, and hold a fixed 
as N goes to infinity. The result is 

C max (a) = -(1 - a) sin (t^-Y (27) 

This equation gives the maximum nearest-neighbor concurrence (under our 
assumptions) for an infinite chain of spin- 1/2 particles in which the overall 
density of up-spins is a. It is reassuring that this formula is identical to the 
one obtained in Ref . || , which considered only infinite chains. Differentiating 
Eq. (pTD , one finds that the optimal value of a is 0.300844, for which C max = 
0.434467. This number is thus our candidate for the maximum nearest- 
neighbor concurrence of an infinite chain of qubits (as in Ref. ||). Note that, 
perhaps surprisingly, for rings of 5 and 8 particles, the maximum values of C 
as given in the above table are smaller than the limiting value for an infinite 
chain. This is no doubt because in these cases one is near the "borderline" 
between two different values of p opt , and neither is particularly good. This 
fact also suggests that the cases N = 5 and N = 8 are the best places to 
look for examples in which the maximum concurrence is not achieved by a 
state satisfying our conditions. 

Indeed, by relaxing condition 2, one can achieve higher entanglement for 
N = 5. The state 

|V> = -^ {sin#[ ITTIII) + ■■•]+ cosfl[ |UTU> + •••]}, (28) 

where the ellipses stand for all translations of the given basis state, has a 
nearest-neighbor concurrence C = 0.468 when 9 = 0.302, which is better than 
the value shown in the above table. We have looked for similar numerical 
improvements for N = 6, 7, 8, 9, and 10, in each case relaxing condition 2 
but preserving condition 1, and we have found none (not even for N — 8). 
Thus it is conceivable that our formula gives the true maximum for certain 
values of N, though it does not do so for all values. In any case, it gives us 
a lower bound on the maximum nearest-neighbor concurrence, which we will 
be able to use in the following Section. 
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To close this section, we write down explicitly the neighboring-pair density 
matrix for our optimal state of the infinite chain. In the form ([3]), the matrix 
elements w and x must both be equal to a, the density of up-spins. This 
is because every up-spin is isolated, so that the probability of the pair state 
|tl) is the same as the probability that the first particle has its spin up, and 
similarly for the probability of |J,|). We already have the value of z, namely, 
half the concurrence; so the density matrix is 



f° 








^ 





0.301 


0.217 








0.217 


0.301 





\o 








0.398 ) 



(29) 

This matrix is not quite one of the special states identified by Munro et al. 
PH , which maximize entanglement for a fixed purity of the density matrix. 
Such a state would have all three of the non-zero diagonal elements equal 
to 1/3. The fact that it is not the same shows that our problem is not 
equivalent to the fixed-purity problem. Nevertheless, it is interesting that 
the two results are as similar as they are. 

3 Entanglement in an Antiferromagnetic Ring 

Though we introduced an effective Hamiltonian in order to solve the preced- 
ing problem, the problem itself did not specify any Hamiltonian. We now 
consider a more concrete physical model of a ring of iV qubits, namely, an 
antiferromagnetic ring of spin- 1/2 particles in which neighboring particles 
interact via the Heisenberg Hamiltonian 

N 

# = £<? 4 -cr m . (30) 

i 

Here a = (a x , a y , <j z ) is the vector of Pauli matrices and, as before, the sum 
i + 1 is understood to wrap around to 1 when i = N. This model has been 
studied extensively over many decades, much of the foundational work having 
been done by Bethe in the early days of quantum mechanics |16|]. In the spirit 



of Section 2 we ask a new question about the model: does the ground state 
maximize the nearest-neighbor entanglement? We restrict our attention to 
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rings with an even number of particles, partly because the calculation is 
considerably simpler in that case, and partly because the symmetry of the 
even-iV ground state suggests an interesting refinement of our question, as 
we will see shortly. 

For the antiferromagnetic ring there is good reason to expect a connection 
between minimizing the energy and maximizing the entanglement. Contrary 
to what one would expect classically, the ground state is not simply the 
alternating state ITiT-l ''•)■ Though this alternating state minimizes the 
energy due to the a z part of the Hamiltonian, it does not do so well for the 
a x and a y parts. By contrast, the ground state for N = 2, which is the singlet 
state 

M = ^(IU>-IIT», (31) 

treats all directions of space equivalently since it is rotationally invariant. 
Intuitively, one expects that for a ring of N particles, each pair of nearest 
neighbors is "trying" to be in the singlet state in order to minimize its own 
energy but is thwarted to some extent by the similar efforts of neighboring 
pairs. Now, the singlet state is maximally entangled; so in a certain sense 
each pair of nearest neighbors, by trying to minimize its energy, is also trying 
to be entangled. We want to see whether the pairs go as far in this direction 
as they possibly could, that is, whether they in fact maximize the nearest- 
neighbor entanglement. Though we do not yet know the maximum possible 
value of this entanglement (because of the extra conditions we imposed on 
our states in Section 2), we can nevertheless use the result of Section 2 as 
a benchmark for evaluating the entanglement of the antiferromagnetic ring. 
For example, if the nearest-neighbor concurrence of the infinite chain is less 
than 0.434467, we know that the entanglement is not maximal. 

We begin by invoking some basic facts about the ground state of an anti- 
ferromagnetic ring with an even number of particles JT7]: it is translationally 



invariant, and it is an eigenstate of the total z-component of spin with eigen- 
value zero. These properties guarantee that the density matrix of each pair 
of neighboring particles has the form 



f v 0\ 

w z 

z w 

V v J 
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(32) 



Let E be the ground-state energy of the system, so that E/N is the con- 
tribution from a single pair: E/N = Tr [p(a ■ a)}. We now re-express the 
energy E/N in terms of the matrix elements of p. The matrix a ■ a, written 
explicitly in the standard basis, is 



a ■ a 



f 1 





\ 





-1 


2 





2 


-1 


\o 





1 j 



(33) 



Thus 

E/N = Ti [p(a- a)] = 2(v - w + 2Rez) = A{v + Rez)-1, (34) 

where we have used the fact that Tr p = 1 . 

It is useful at this point to write the matrix elements v and z in terms of 
the coefficients that define the ring's state \ip). Just as in Section 2, we can 
write \ip) as 

\^)= Y. b h ...i p \ii---i P ), (35) 

l<il<-<i p <JV 

where p now has the specific value N/2. And just as before, we have 

z = z2 biM-kph+iM-kp- (36) 

l<k 2 <-<k p <N 

The corresponding expression for the matrix element v is 

v = H |6»,i+i,fc3...*p| 2 - (37) 

l<fc 3 <-<fc p <A r 

Note that changes in the phases of the coefficients b do not affect v, though 
they do affect z. In order to minimize the energy as given in Eq. (|3~4]), we 
want to choose these phases so that Rez is as negative as possible. For a 
fixed set of absolute values of the 6's, this can be done be letting all the fe's 
be real, with alternating signs given by 

sign of b n .., p = (-l)^+-+v (38) 

in which case every term of Eq. (|36|) is negative or zero. Thus for the ground 
state of this system, we can write the energy per particle as 

E/N = A{v-\z\)-l. (39) 
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Now, recall that the concurrence of a density matrix of the form (|32| ) is 

[Eq- (1)] 

C = max {2(|z| -v),0}. (40) 

We thus arrive at the following expression for the concurrence C gs of the 
ground state of this system, assuming (as is the case) that the ground-state 
energy is sufficiently negative to make C gs positive. 



a 



.</* 



-\(E/N) + 1}. 



(41) 



This simple relationship depends on the fact that the number of particles in 
the ring is even. If N were odd, the pair density matrix would not have the 
form (^) and its concurrence would most likely not be a simple function of 
the energy alone. 

The ground state energies of antiferromagnetic rings have been computed 
for many values of N |L7, 18 1, including the limiting case N 



OC 



191. From 



these results and Eq. ( |4"T| ) we can immediately write down the concurrences. 
The following table shows the values of C gs for several values of N, along 
with corresponding values of C max that we computed in Section 2. The 
figure 0.386 appearing in the table as the concurrence of the ground state of 
the infinite chain can be written exactly as 2 In 2 — 1. 



N 


2 


4 


6 


8 


10 


oo 


-E/N 


3.000 


2.000 


1.868 


1.825 


1.806 • 


• 1.773 


Cg s 


1.000 


0.500 


0.434 


0.412 


0.403 • 


• 0.386 


KJ max 


1.000 


0.500 


0.471 


0.433 


0.449 • 


■ 0.434 



Thus, though for very small rings the antiferromagnetic ground states are as 
entangled as the states we found in Section 2, for larger rings they fall short. 
We can therefore conclude that the ground state of an antiferromagnetic ring 
does not in general maximize nearest-neighbor entanglement. 

There is, however, a more limited sense in which these ground states do 
maximize entanglement; this is the refinement we mentioned earlier. Let us 
restrict our attention to the set of states which, like the ground state, are 
translationally invariant and have zero total z-component of spin. We will 
call such states "balanced." We now show that the antiferromagnetic ground 
state maximizes entanglement within the set of balanced states. 
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Let us divide the set of all balanced states into equivalence classes, two 
states being called equivalent if their coefficients bi x i in Eq. (|35|) agree 
in magnitude, differing only in their phases. Of all the states in a given 
equivalence class, none has a greater nearest-neighbor concurrence than the 
unique state in that class for which the phases are given by Eq. fl38|) . This 
is because, as in the case of the ground state, such phases allow perfect 
constructive interference in Eq. (^). To put it in symbols, C(ip) < C(i/j ), 
where \ip) is a general balanced state and \ipo) is the state obtained from 
\ip) by adjusting the phases of the 6's in accordance with Eq. (|38|). Now, 
for |^o), the expectation value of the energy (ip \H\ipo) is given by the same 
expression as in Eq. (^): 



-(1>o\H\1h>)=4(v-\z\)-1. 



(42) 



The concurrence of |^o) is given by Eq. ([TOP; so we have 

W) < C{,/> ) = max {-| [i^olWo) + l] ,0J 
< max{-±[(E/N) + l],0} = C g8 . 



(43) 



The last inequality comes from the fact that the ground state minimizes the 
expectation value of the energy. We have thus shown that no balanced state 
has a nearest-neighbor concurrence larger than that of the ground state. 

For comparison with Eq. (^9|) , it is interesting to write down explicitly the 
neighboring-pair density matrix for the ground state of an infinite antiferro- 
magnetic chain. This density matrix is uniquely determined by the value of 
the concurrence, C = 2 In 2 — 1, and by the fact that the state is rotationally 
invariant (the latter condition implies that \z\ + v = w). One finds that 



(44) 



If we think of the spins of the antiferromagnetic chain as "trying" to max- 
imize their entanglement, then evidently they are using a rather different 
strategy than the one we used in Section 2. There is no longer any pro- 
hibition against neighboring up-spins. Indeed the presence of such up-spin 



f 0.102 








° ^ 





0.398 


-0.295 








-0.295 


0.398 





V o 








0.102 / 
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pairs in the antiferromagnetic chain allows the off-diagonal element z to have 
a larger magnitude than in Eq. (p9j ), which is good for entanglement. On 
the other hand, the presence of such pairs also forces the matrix element 
v = (TTIpITT) to be non-zero, which is what reduces the concurrence to a 
value less than our best value of Section 2. 



4 Conclusions 

We have obtained two main results. 

First, for a ring of iV qubits in a translationally invariant state, we have 
found values of the nearest-neighbor concurrences that we know to be achiev- 
able and that for some values of N may be optimal. At the least, they are 
lower bounds on the maximum possible concurrences. 

Second, we have found that the ground state of an antiferromagnetic ring 
with an even number of particles typically does not maximize the nearest- 
neighbor concurrence over all states, but that it does achieve such a maximum 
over the set of translationally invariant states with no net spin in the z 
direction. This set of "balanced" states includes all the eigenstates of total 
spin with eigenvalue zero; so we can also say that the ground state maximizes 
C relative to all the spin-0, or rotationally invariant states. 

Putting Sections 2 and 3 together, we can conclude that whatever the 
maximum-concurrence states may be, they are certainly not balanced. In 
other words, for maximizing concurrence it is best to have one direction of 
spin favored over the opposite direction. This is perhaps counterintuitive, 
since a maximally unbalanced state such as |'j 1 1 • • •} is not entangled at all. 

The subject of Section 3 represents an unusual mix: one does not often as- 
sociate entanglement with energy-minimization. One might wonder whether 
the entanglement-maximization property of antiferromagnetic rings, limited 
though it is, is a special case of a more general connection between energy 
and entanglement. Do physical systems tend to favor entangled states over 
unentangled states? In a straightforward interpretation of this question, the 
answer would seem to be no. Ferromagnetic systems, for example, have 
ground states in which the spins are completely unentangled. Perhaps one 
could identify a special class of Hamiltonians with interesting entanglement- 
maximizing properties, but at present it is not clear how large such a class 
might be. 



17 



We would like to thank Daniel Aalberts for a number of very helpful 
discussions. 



References 

[1] See, for example, C. H. Bennett, H. J. Bernstein, S. Popescu, and 
B. Schumacher, Phys. Rev. A 53, 2046 (1996); V. Vedral, M. B. Plenio, 
M. A. Rippin, and P. L. Knight, Phys. Rev. Lett. 78, 2275 (1997); M. 
Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 

[2] B. Schumacher, Phys. Rev. A 54, 2614 (1996). 

[3] D. P. DiVincenzo, Science 270, 255 (1995). 

[4] J. S. Bell, Physics 1, 195 (1964). 

[5] V. Coff man, J. Kundu, an d W. K. Wootters, Phys. Rev. A 61, 052306 
(2000); huant-ph/9907047| . See also D. BruB, Phys. Rev. A 60, 4344 



(1999). 
[6] W. Diir, |quant-ph/0006T05 



[7] M. Koashi, V. Buzek, and N. Imoto, |quant-ph/0007086 . 

[8] W. K. Wootters, huant-ph/OOOTTTl 

[9] See, for example, A. V. Thapliyal, Phys. Rev. A 59, 3336 (1999); J. 
Kempe, Phys. Rev. A 60, 910 (1999); D. Aharonov, |quant-ph/9910081 



C. H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin and A. V. Thap- 
liyal, |quant-ph/9908073| ; N. Linden, S. Popescu, B. Schumacher and M. 



Westmoreland, |quant-ph/9912039| ; G. Vidal, W. Diir and J. I. Cirac, 



Phys. Rev. Lett. 85, 658 (2000). 

[10] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 

[11] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[12] C. H. Bennett, D. P. DiVincenzo, J. Smolin, and W. K. Wootters, 
Phys. Rev. A 54, 3824 (1996). 

[13] S. Ishizaka and T. Hiroshima, Phys. Rev. A 62, 22310 (2000). 

[14] W. J. Munro, D. F. V. James, A. G. White, and P. G. Kwiat, unpub- 
lished. 



19 



[15] E. Lieb, T. Schultz, and D. Mattis, Annals of Physics 16, 407 (1961). 

[16] H. A. Bethe, Z. Physik 71, 205 (1931). 

[17] R. L. Orbach, Phys. Rev. 115, 1181 (1959). 

[18] J. C. Bonner and M. E. Fisher, Phys. Rev. 135, A640 (1964). 

[19] L. Hulthen, Arkiv. Mat. Astron. Fysik 26A, No. 1 (1938). 



20 



